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The rich internal structure and long-range 
dipole-dipole interactions establish polar 
molecules as unique instruments for quantum- 
controlled applications and fundamental in- 
vestigations. Their potential fully unfolds at 
ultracold temperatures, where a plethora of 
effects is predicted in many-body physics [!Q [2], 
quantum information science p£l [4], ultracold 
chemistry [5, 6j, and physics beyond the standard 
model [7, 8j. These objectives have inspired the 
development of a wide range of methods to 
produce cold molecular ensembles [9l-fl4j . How- 
ever, cooling polyatomic molecules to ultracold 
temperatures has until now seemed intractable. 
Here we report on the experimental realization of 
opto-electrical cooling [15j, a paradigm-changing 
cooling and accumulation method for polar 
molecules. Its key attribute is the removal of 
a large fraction of a molecule's kinetic energy 
in each step of the cooling cycle via a Sisyphus 
effect, allowing cooling with only few dissipative 
decay processes. We demonstrate its potential by 
reducing the temperature of about 10 6 trapped 
CH 3 F molecules by a factor of 13.5, with the 
phase-space density increased by a factor of 29 
or a factor of 70 discounting trap losses. In con- 
trast to other cooling mechanisms, our scheme 
proceeds in a trap, cools in all three dimensions, 
and works for a large variety of polar molecules. 
With no fundamental temperature limit antic- 
ipated down to the photon-recoil temperature 
in the nanokelvin range, our method eliminates 
the primary hurdle in producing ultracold poly- 
atomic molecules. The low temperatures, large 
molecule numbers and long trapping times up to 
27 s will allow an interaction-dominated regime 
to be attained, enabling collision studies and 
investigation of evaporative cooling toward a 
BEC of polyatomic molecules. 

The ability to prepare ultracold molecular ensembles 
has an application potential akin to that of ultracold 
atoms some decades ago. In fact, the association of 
KRb dimers [16 as well as the laser cooling of SrF [17] 
has brought fascinating physics within reach. However, 
both approaches are restricted to a highly specialized set 
of purely diatomic molecule species. In order to inves- 
tigate fundamental physics based on relativistic effects 
near heavy nuclei or parity violation effects in chiral 
molecules, or to study molecules of astrophysical, bio- 
logical, or chemical interest, a more general approach 



to preparing ultracold molecular ensemble is imperative. 
This holds in particular for the rich chemical variety of 
carbon-, nitrogen-, or oxygen-based molecules for which 
the constituent atoms have not even been laser cooled. 
Devising a dissipative process to cool such molecules into 
the ultracold regime has been an exceedingly challeng- 
ing problem. The standard approach for atoms, laser 
cooling, is in general impossible for molecules due to the 
lack of suitable cycling transitions. Creating an artificial 
cycling transition via cavity cooling [18] has not been 
demonstrated despite substantial experimental [19] and 
theoretical [20-22J effort. Likewise, evaporative or sym- 
pathetic cooling to ultracold temperatures [23] has not 
been realized due to lack of density or losses from inelas- 
tic collisions. 

A particularly promising general framework to cool 
molecular ensembles is to replace the weak photon re- 
coil in laser cooling with sufficiently strong forces to re- 
move the entire molecule's kinetic energy in a single pro- 
cess [T5l [24] . A first step towards this approach is the re- 
cent achievement of accumulation of NH molecules using 
a single-photon transition [25 . Here we present a full im- 
plementation of opto-electrical cooling [15 for molecules 
stored in an electric trap, featuring both accumulation 
and cooling. Energy is extracted by allowing molecules 
to move up and down an electric field gradient in different 
states with differing Stark energies while dissipation to 
remove entropy is provided by a spontaneous vibrational 
decay. The scheme is applicable to all molecules with 
strong electric field interaction and pure ro-vibrational 
states, and thus constitutes a general method for cooling 
molecules to ultracold temperatures. 

Reducing the temperature by a large factor requires 
the cooling cycle to be repeated, and consequently, con- 
trol being maintained over the internal molecular state. 
This is achieved using the level scheme shown in Fig. [ljt. 
An electric trap as depicted in Fig. [TJ) with a homo- 
geneous electric field in the trap centre and strongly 
increasing fields near the trap boundary [26 leads to 
a position-dependent potential energy for the relevant 
molecular states as shown. We label states with vi- 
brational quantum number v and symmetric-top rota- 
tional quantum numbers J, M [27 as \v\ J, ^K, ±M) 
with =fK chosen positive. For a parallel vibrational 
transition with AK = 0, the state |1;3,3,3) decays to 
the four rotational v = states shown, with decay to 
the |0;4,3,2) state ignored due to the small Clebsch- 
Gordan coefficient of Coupling the more weakly 

trapped |0;3,3,2) and |0;4, 3,3) states with microwave 
(MW) radiation and driving the |0;3,3,2) to |1;3,3,3) 
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FIG. 1: Implementation of molecule cooling, a, Radiation couples molecular states which experience a position-dependent 
potential energy in an electric trap as indicated. Molecules following the emphasized route lose more kinetic energy when 
entering the strong-field edge region of the trap than they regain when returning to the trap centre in a more weakly trapped 
state. Combined with the unidirectionality of the optical pumping back to the strongly trapped states, this leads to cooling, 
b, During cooling the molecules are confined in a microstructured electric trap [26] , with the IR, MW, and RF fields applied 
as indicated, c, Experimental sequence for cooling: warm molecules are continuously loaded into the trap for 16 s with the IR, 
MW, and two RF fields at 3.4 GHz and 1.6 GHz already applied. The trap is then electrically closed [26] and the RF frequency is 
decreased stepwise to 390 MHz over the next 21 s. Before unloading the molecules for detection, an internal-state-discrimination 
(ISD) technique is applied during 4 additional seconds (see methods). 



transition with an infrared (IR) laser results in optical 
pumping to the strongly trapped |0; 3, 3, 3) and |0; 4, 3, 4) 
states. Adding a radio-frequency (RF) to couple neigh- 
bouring M-sublevels in strong electric fields completes 
the opto-electrical cooling cycle. Losses to the untrapped 
M = states are avoided by coupling the neighbouring 
M-sublevels with the RF at a rate which is slow com- 
pared to the optical pumping [28 j. Note that Stark de- 
tuning due to the electric fields plays a key role in se- 
lectively addressing only the desired IR and MW transi- 
tions [28] . 

The experimental sequence to prepare and detect a 
sample of cooled molecules is shown and explained in 
Fig. In addition to cooling, we apply the following 
two experimental sequences. First, a reference no-cooling 
measurement with the IR, MW, and RF fields left off and 
the storage time shortened by 20 s results in an uncooled 
molecular ensemble. Second, reintroduction of the IR 
and MW fields results in accumulation of molecules, also 
without cooling. Here, molecules entering the trap in the 
weakly trapped states |0;3,3,M) and |0;4,3,M) with 
1 < M < J — 1 are pumped to the strongly trapped 



|0;3,3,3) and |0;4,3,4) states. 

All experimental sequences end by unloading the 
molecules from the trap and guiding them to the ioniza- 
tion volume of a quadrupole mass spectrometer (QMS) 
for detection. This results in signals as shown in the in- 
set to Fig. [TJ), in this case for the full cooling cycle to 
390 MHz as well as for the same cycle but with the IR 
laser left off. Remarkably, the sequence with laser on re- 
sults in a strong molecule signal whereas leaving the laser 
off eliminates the signal almost entirely. Thus, the exci- 
tation by the laser and the subsequent spontaneous decay 
clearly plays a key role in keeping molecules in the trap, 
indicating that a key element of the cooling cycle is work- 
ing. Note that state-discriminating measurements (see 
methods) show that 80 % of the molecules with laser on 
are in the states |0; 3, 3, M) and |0; 4, 3, M) with most of 
these in the top M-sublevels [28] , establishing our exper- 
iment as a practically state-selected source of molecules. 
This compares to about 20% of molecules initially en- 
tering the trap in the states |0;3,3,M) and |0;4, 3, M) 
potentially involved in the cooling. 

For proof of cooling, we now present four complemen- 
tary measurements. First, we examine the effect of the 
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FIG. 2: Demonstration of opto-electrical cooling, a, Detected molecule signal vs. trap electric field strength during 
unloading. In addition to results for the no-cooling measurement, the accumulation, and the full cooling cycle to 390 MHz, 
we show results for cooling to the intermediate frequency of 1.22 GHz. Cooling results in a more than 10-fold decrease in 
the optimal unloading voltage and, together with the accumulation, in a multifold signal increase at low unloading voltages, 
b, Relative unloading signal with RF as a knife-edge filter applied. The legend shows the RF cooling frequency (except for 
accumulation) after which the much stronger RF knife is applied. We subtract 180 MHz from the RF knife frequency to account 
for the offset electric field at the centre of the trap. 



trap electric field strength on the molecules during un- 
loading. The field strength during this time is partic- 
ularly important since it determines both the confining 
potential for the molecules and the molecule transport 
efficiency to the detection QMS. In particular, a lower 
electric field strength reduces the trap depth so that hot- 
ter molecules are lost, and a higher electric field strength 
during unloading inhibits colder molecules from finding 
the trap exit hole to the QMS. We thus expect that 
the molecule signal falls off in both directions from a 
maximum at a confining electric field strength which is 
proportional to the temperature of the molecules, de- 
pending on the molecule Stark shifts. Fig. |2^l shows 
the measured molecule signal vs. unloading electric field 
strength for four different experimental cycles. Applica- 
tion of the cooling sequence results in a strong decrease in 
the optimal unloading electric field strength, from about 
45kV/cm to 16kV/cm and 3.3kV/cm. Since the plotted 
signal is mainly from molecules in the states |0;3,3,3) 
and |0;4, 3,4) [28] which have a similar Stark shift, our 
data clearly shows that a large amount of cooling takes 
place. The relatively low signal for the full cooling cy- 
cle is to a large part due to trap losses during the ad- 
ditional 20 s of storage compared to accumulation and 
no-cooling, accounting for a 2.4 times smaller signal (see 
below). Note that the accumulation alone results in ap- 
proximately a 2-fold larger molecules signal, thus proving 
an increase in the phase-space density already here. 

As a second measurement, we use a strong RF field 
as a knife to drive transitions to untrapped states and 
thus eliminate hotter molecules from the trap. Results 



are shown in Fig. [2J3- For all three measurements with 
cooling, the cooling cycle is completed to 390 MHz. How- 
ever, two additional seconds of storage are inserted im- 
mediately after the indicated cooling frequency for ap- 
plication of the RF knife. This allows for practically 
identical experimental sequences. Only for accumula- 
tion, the molecules are unloaded directly. We observe 
a sharp flank in the molecule signal at a knife frequency 
which strongly depends on the frequency of the cooling 
step just before application of the knife. The position of 
the flank is shifted by approximately a factor of 10, sim- 
ilar to the shift observed in Fig. [2^l. Again, the dramatic 
shift can only be explained by a large amount of cooling. 
Note that the flanks of the four measurements are well 
separated, indicating that the energy distributions of the 
molecules in the four measurements hardly overlap. 

Both of the previous measurements are sensitive to the 
internal molecular state since Stark shift and Stark split- 
ting directly determine the effect of the trap voltages 
and RF frequencies, respectively. To avoid this effect, 
we perform a third experiment and directly measure the 
velocities of the molecules via a time-of-flight technique. 
Specifically, we switch the 55 cm long quadrupole guide 
segment leading to the QMS on and off and measure the 
arrival time of the molecules. Note that the unloading 
voltages affect the distribution of molecules arriving at 
the QMS, with colder (hotter) molecules being preferen- 
tially detected at lower (higher) voltages [28]. However, 
by measuring at the optimal unloading voltages, the ef- 
fect will be the same for all measurements relative to the 
ensemble in the trap. We use the same four experimental 
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FIG. 3: Time-of-flight measurements to determine the 
molecule velocities, a, Rising edge of the normalized trap 
unloading signal for the same variations of the experimental 
sequence as in Fig. [2Jl For cooled molecules one needs to 
wait over three times longer before molecules arrive at the 
QMS. b, Molecular velocity distributions derived from the 
rising-edge signals shown in (a) taking into account the num- 
ber of detected molecules. The dashed blue curve shows the 
hypothetical result if trap losses (shown in Fig. [4| could be 
eliminated. The smooth appearance of the data is due to con- 
volution with a Gaussian. The shaded regions represent the 
la statistical error. We include the mean velocity and tem- 
perature calculated from the data sets, c, Measured molecule 
temperature and total change in phase-space density (PSD) 
after each intermediate RF frequency as well as for no-cooling 
and accumulation. 
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FIG. 4: Trap lifetime for cooled and uncooled 
molecules. Cooling substantially increases the trap lifetime. 
The inset shows the raw data for three of the lifetime mea- 
surements. 



sequences as for Fig. |2^i with the results shown in Fig. [3] 
Calculating the mean longitudinal velocity (v z ) from 
the data and deriving the molecular temperature T via 
\ksT = \m(v z ) 2 we arrive at the central result: the 
molecules have been cooled by a factor 13.5 from 390 mK 
to 29 mK. Together with the estimated molecule density, 
this corresponds to a factor 29 increase in phase-space 
density. In fact, the cooling and accumulation by them- 
selves increase the phase-space density by a factor of 70, 
but this is reduced by trap losses during cooling. 

As a fourth and final measurement, Fig. [4] shows the 
lifetime of molecules in the trap as a function of the final 
cooling frequency. As expected for cooling [26 , the life- 
time increases from about 10 s without cooling to about 
27 s with maximal cooling. This dramatic change con- 
firms once more that cooling takes place. By integrating 
the measured loss rate over time, we find that the addi- 
tional 20 s of storage for cooling to 390 MHz compared to 
no-cooling or accumulation leads to a loss of about 60 % 
of the molecules, as has been used above in the discussion 
of Fig. [2^i and in Fig. (3}d. For better trapping or faster 
cooling, this loss could be reduced. 

In principle, opto-electrical cooling can be continued to 
much lower temperatures. However, losses during cool- 
ing and the low detection efficiency make further cool- 
ing increasingly tedious for the present setup. Here, 
an improved trap design could increase the trap life- 
time, and higher detection efficiencies would make losses 
more tolerable. Varying the cooling cycle or choosing a 
molecule species with a faster spontaneous decay would 
allow faster cooling. 

Even without these improvements, the present exper- 
iment brings key applications of cold polar molecules 
within reach. First, the increased phase-space density 
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improves the sensitivity for high-resolution spectroscopy 
and collision experiments. Second, the obtained tem- 
perature allows loading of molecules into more weakly 
confining microwave [29] or optical traps which can hold 
molecules in their rotational ground state, a prerequisite 
for achieving quantum degeneracy. Finally, our chip-like 
trap and guide architecture matches the demands on hy- 
brid quantum systems for quantum information process- 
ing with cold and ultracold molecules [4] . 

Methods 

The initial sample of molecules is generated by ve- 
locity filtering [12] from a liquid-nitrogen cooled source 
and loaded into the electric trap [26 via a quadrupole 
guide. Detection is performed by guiding the molecules 
to the QMS via a second guide. The selection of CH3F 
as molecule species is based on favourable properties for 
trap loading and detection. Large rotational constants 
and sufficient vapour pressure down to almost 100 K al- 
low efficient velocity filtering with an adequate fraction 
(~ 20 %) of molecules in the rotational states used for 
cooling. Low contributions from other molecules at the 
atomic mass 34 of CH3F results in very low background 
for the QMS detection. For dissipation, the approxi- 
mately 15 Hz spontaneous decay rate of the most suited 
parallel vibrational v\ symmetric C-H stretch mode is 
sufficient. A description of the IR and MW radiation 



sources required for CH 3 F [30] is provided in the ap- 
pendix. 

For internal- state-discriminating measurement, a mi- 
crowave depletion (MWD) pulse consisting of various 
MW frequencies near 204 GHz is applied during the last 
four seconds before trap unloading. This mixes all M- 
sublevels of the J = 3 and J = 4 states, leaving only the 
unaffected background from molecules in other states in 
the trap. The difference in the unloading signal with and 
without the MWD applied is used for all measurements 
shown in Figs. [2}|4] 

The offset electric field in the centre region of the trap 
results in an offset Stark splitting of 180 MHz between 
neighbouring M-sublevels for molecules with J = \K\ = 
3. For molecules with J = 4, all Stark splittings are 
smaller, but this is inconsequential due to microwaves 
being used to couple the |0; 3, 3, M) and |0; 4, 3, M + 1) 
states: the RF transitions will predominantly occur for 
J = 3. This is particularly important for the RF filter 
measurement in Fig. [2}d. 
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Appendix A: Detailed Cooling Scheme 



The illustration of the opto-electrical cooling cycle in Fig. [TJi in the main text includes only a simplified level scheme 
consisting of the molecular states and transitions which constitute the most elementary cooling cycle. Fig. [5] shows 
an expanded version of the cooling cycle where we consider additional M-sublevels of the J = 3 and J = 4 states as 
well as additional decay channels from the v = 1 state. We find that including the lower trapped M-sublevels in the 
cooling cycle leads to a substantial reduction in losses to untrapped M = states and that the additional spontaneous 
decay channels for the most part do not affect the cooling cycle. 
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FIG. 5: Detailed energy-level diagram for cooling showing the complete set of IR, MW and RF transitions as well as 
all spontaneous decay channels involved in cooling. Black dashed arrows indicate the key spontaneous decay channels which 
transfer the molecules back towards the maximally trapped states, with all other decay channels indicated by gray dashed 
arrows. The dashed energy levels are loss channels. 

Due to the constant differential Stark splitting between neighbouring M-sublevels for J and K fixed, it is basically 
impossible to drive RF transitions from a state M = J to M = J — 1 (as required for cooling) without also driving 
transitions to lower M-sublevels and ultimately to untrapped M = states. For the cooling cycle, these losses are 
reduced by driving the RF transitions at a slow rate compared to the optical pumping back to the M = J states. 
For an additional substantial reduction of these losses, the states |0;3,3, 1) and |0;4, 3,2) can be incorporated into 
the cooling cycle. Coupling these states via microwaves and coupling to the |1; 3, 3, 2) state via the IR laser results in 
optical pumping back to the core states of the elementary cooling cycle, as shown in Fig. [5j This doubles the number 
of "unlucky" AM = — 1 RF transitions which need to occur before a molecule is lost. 

Fig.[]Ji in the main text only shows spontaneous vibrational decay to the |0; 3, 3, 3) and |0; 4, 3, 4) states. Obviously, 
spontaneous decay to other states will also occur. However, since spontaneous vibrational decay with A J = 0, AM = 
— 1 and A J = +1,AM = populates states which are strongly coupled to the v = 1 state via the IR and MW, 
the effect on the molecules is negligible. This leaves the decay with A J = +1, AM = —1 as the only decay channel 
which can induce losses. However, the Clebsch-Gordan coefficient of the |1;3,3,3) to |0;4, 3,2) and the |1;3,3,2) to 
|0; 4, 3, 1) transition is only and respectively. Moreover, molecules decaying to the |0; 4, 3, 2) state will generally 
be repumped to higher M states via the |1; 3, 3, 2) state. Considering the few number of spontaneous decays required 
for cooling, losses via the A J = +1, AM = —1 spontaneous decay channel can therefore be ignored for the present 
experiment. 

As a final consideration, we discuss two possibilities for spontaneous vibrational decay to states other than |0; 3, 3, M) 
and |0;4, 3,M). First, due to the electric fields in the trap, molecular states with different J are coupled so that 
spontaneous decay with A J > 1 in principle becomes possible. However, this was analyzed for the original proposal of 
opto-electrical cooling [15] where it was found that even for CF3H with much smaller rotational constants, spontaneous 
decay with A J = 2 is marginal even for fields as large as lOOkV/cm. Second, it is well known that the v\ C-H 
vibrational stretch mode of CH3F is coupled to the 2v2 and 2v$ modes via a strong Fermi-resonance [3TJ [32]. As 
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a result, spontaneous decay to the first excitation of the ^5 mode will occur, with a rate of about 1 — 2 Hz [33] . 
Thus, approximately every 10th decay will occur via this channel. The Fermi resonance thus causes a significant loss 
channel, evidence for which is discussed in appendix [D] 



Appendix B: Rate estimates for opto-electrical cooling 

In this section we perform some simple back-of-the-envelope type estimates for the dynamics of opto-electrical 
cooling. These estimates allow us to derive a scaling law for the RF losses to M = as well as to explain the 
discrepancy between the vibrational spontaneous decay rate of about 15 Hz and the cooling time of about 8 s per 
halving of the molecular temperature. 

Our entire analysis of opto-electrical cooling relies on using rates to describe transitions between populations in 
various molecular states. We thereby ignore the possibility of any coherences, an assumption which is strongly justified 
by the experimental parameters. In particular, the RF transition rate of about 1 Hz during cooling compares to a 
trap frequency of about 1 kHz: the molecules thus come into resonance with the RF many times before a transition 
occurs. Similarly, the MW and IR rate of about 100 Hz during cooling compares to a Stark broadening of many MHz 
in the central region of the trap. 

For the following discussion, it is useful to group the molecules into three sets of states, the strongly trapped pair of 
states |0;3,3,3) and |0;4,3,4), the (on average) moderately trapped set of states |0;3,3,2), |0;4, 3,3), and |1;3,3,3), 
and the (on average) weakly trapped set of states |0;3,3, 1), |0;4, 3,2), and |1;3,3,2). These sets are motivated by 
the fact that the molecular states in each of these sets are strongly coupled to each other by IR and MW radiation. 
As a result, we can treat each set as a single doubly- or triply-degenerate state. 

Based on the previous grouping of molecular states, cooling can be modelled as an RF which couples neighbouring 
sets of states with a rate of 7rf, and a repumping to more strongly trapped states with a rate of 71R. With 
populations p s , p m , and p w in the strongly, moderately, and weakly trapped set of states, respectively, and the 
assumption 7rf <C 7ir, we find an approximate steady state for cooling with 7rfPs = JmPm and 7rf Pm — 7irPw 

We therefore find that p w = p s x • Since 7rf is simply the rate of cooling and with losses to untrapped states 

given by 7rfPw, we find that the loss rate to untrapped states is suppressed compared to the cooling rate by a factor 

of ( ) . Note that without the weakly trapped set of states, the second factor of would be absent and losses 

\ 7IR / 7IR 

would scale much less favourably for small 7rf- Also, we have found that to avoid losses, the assumption 7rf <C 7ir is 
true self consistently. This is equivalent to p w <C p m <C p s , and we have therefore found that the majority of molecules 
must be in the strongly trapped states at any given time. 

We now estimate realistic rates for our experiment. The spontaneous decay rate of the vibrationally excited state is 
approximately 15 Hz. However, the population in both the moderately and weakly trapped set of states is spread over 
three states, only one of which is in v = 1, effectively reducing 71R by a factor three. Moreover, the branching ratio 
of decay from the state |1;3,3,3) to the strongly trapped states is about 0.75, further reducing 71R by this factor. 
We therefore find 71R = 4 Hz, and to sufficiently suppress losses, applying the RF with a rate of 7rf = 1 Hz seems 
reasonable. 

To obtain the rate at which the temperature can be reduced, two additional effects play a role. First, we estimate 
the number of cooling cycles required to achieve a given reduction in temperature. Consider a molecule with energy 
E in the rotational sublevel M. The condition for the molecule to have sufficient energy to reach resonance with an 
RF frequency of /rf (allowing an RF transition to take place) is E > M/rf- Since the RF transition reduces the 
molecule's energy by a factor of E/(E — /rf), optimal cooling takes place for /rf as large as possible. To account 
for a range of molecule energies, we consider /rf such that one final RF transition is possible. For M = 3, this is 
equivalent to 3 /rf < E < 4 /rf which translates to an average reduction in energy by a factor of about 3.5/2.5 = y/2. 
We thus find that approximately two cooling cycles are necessary to reduce the temperature by a factor of 2. 

The final effect for the cooling duration is that due to the stochastic nature of the cooling process, a cooling rate 
of 1 Hz means that a fraction of e _1 of the molecules will remain uncooled after Is. The cooling rate 7rf is thus 
substantially faster than the time required to efficiently cool molecules by one step. For a reasonable cooling efficiency 
of 1 — e -3 , this adds a further factor of 3 to the cooling time. We therefore end up with a total of 6 s required to 
cool the molecules by a factor of 2, which, within the accuracy of our estimate, agrees with the 8 s chosen for the 
experiment. 
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FIG. 6: Accumulation. We show the total integrated QMS signal during unloading as a function of the time during which the 
radiation fields for accumulation are applied. The solid line is the expected time dependence of the signal from the rate-equation 
model, showing excellent agreement with the experimental data. 



Appendix C: Rate equations 

In the previous section we showed that even a very simple estimate based on rates allows several aspects of opto- 
electrical cooling to be understood quantitatively. In this section, we expand our model by considering populations of 
molecules in the trapped states shown in Fig. [5] and solving time-dependent rate equations describing the transitions 
between the various states. Due to the ability to independently determine all relevant parameters other than the 
absolute signal amplitude, the expanded model works particularly well for an "accumulation" measurement where 
we perform optical pumping from the moderately to the strongly trapped states. The perfect fit of the experimental 
data with the values expected from our model demonstrates that our model accurately describes the processes in the 
trap and, moreover, that accumulation works as intended. For cooling, determining all the relevant parameters for 
the model is not possible, resulting in only an approximate description of our data. However, the model nonetheless 
provides additional insight into the cooling process. 

For the rate-equation model, we consider the time-dependent populations of molecules Pi(t) with i = 1, ...,9 in the 
seven low-field seeking states |0;3,3,M) and |0;4, 3, M) as well as in the two excited states |1;3,3,3) and |1;3,3,2). 
The state |1;3,3, 1) can be ignored since it is strongly coupled to the untrapped |0;3,3,0) state via the IR laser. 
For a given initial state distribution Pi(0), we determine the state evolution by solving the set of rate equations 
Pi(t) = AijPj(t). Here, Aij is the rate at which molecules in the state j are transferred to the state z, describing 
the effect of the radiation fields which are applied during the experiment as well as the spontaneous decay from the 
v = 1 states. 



1. Accumulation 

Our rate-equation model can be applied particularly successfully to an accumulation measurement where the 
molecules are optically pumped from the moderately trapped states |0;3,3,2) and |0;4, 3,3) to the strongly trapped 
states |0; 3, 3, 3) and |0; 4, 3, 4). Experimentally, we start by performing the cooling sequence to 1.22 GHz as described 
in the main text. This is followed by two seconds during which microwaves are applied to couple the four states 
|0; 4, 3, 4), |0; 3, 3, 3), |0; 4, 3, 3), and |0; 3, 3, 2), thereby distributing the molecules among these four states. During the 
subsequent four seconds, we apply the MW and IR fields shown in Fig. [5] (but no RF), identical to the accumulation 
as described in the main paper, for a variable amount of time T. This causes a certain fraction of the molecules to 
be repumped to the strongly trapped states. After the four seconds, we reduce the confining electric field strength of 
the trap and unload the molecules from the trap for detection, with the molecular state affecting the probability for 
detection. 

The accumulation measurement is particularly suited for comparison to the rate-equation model since all aspects 
of the experiment not included in the model can be accounted for. In particular, the energy distribution of the 
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FIG. 7: Cooling efficiency vs. RF power. We show the number of detected molecules as a function of the power of the 
390 MHz RF during the last 4 s of cooling. The solid line is the expected dependence of the signal from the rate-equation model. 



molecules in the trap is not included in the model. The energy of a molecule strongly affects its detection efficiency 
and can in principle affect the transition rates between the molecular states. However, due to the box-like shape 
of our trap potential [26] . the probability for a molecule to be in the central trap region is close to unity and, in 
particular, rather insensitive to a molecule's energy Hence, the rate of the MW and IR transitions, which occur in 
the central trap region, will also be quite insensitive to a molecule's energy. As a result, we expect the two seconds of 
mixing to create similar energy distributions for the molecules in each of the molecular states and for this to remain 
the case throughout the accumulation. Moreover, we expect the energy distribution in each of the final states to be 
independent of the duration of the accumulation. Consequently, we can define energy- averaged detection probabilities 
d s and d m for molecules in the strongly and moderately trapped states, respectively, which are independent of the 
accumulation time T. 

Based on the previous discussion, the detected signal for the accumulation measurement is composed of individual 
contributions as follows. Molecules which are in the strongly trapped states at the start of the accumulation are 
unaffected by the accumulation and, together with the background count rate, will provide a constant contribution 
S s to the detected signal. For molecules initially in the moderately trapped states, a fraction p(T) depending on the 
duration of the accumulation ends up in the strongly trapped states and the rest remains in the moderately trapped 
states. For N m molecules initially in the moderately trapped states, this results in a contribution N m d s p(T) + 
iV m d m (l — p(T)) to the detected signal. As a result, the total signal is of the form a + bp(T), where a and b are 
taken as fit parameters. Since the MW and IR transition rates as well as the spontaneous decay rate are known from 
other measurements, the function p(T) can be precisely calculated, so that the vertical scaling and vertical position 
remarkably remain the only fit parameters. 

The experimental data and the expected time dependence obtained from the rate-equation model is shown in Fig. [6] 
The value of chi-squared for the variation of the data around our model is x 2 = 7.4, which, for our v = 12 degrees 
of freedom results in x 2 /v = 0.617, showing that the model and our data are in perfect agreement. Note that after 
about 1 s, the signal no longer increases. This is due to all molecules having been pumped to the strongly trapped 
states by this time. While accumulation from the weakly trapped to the strongly trapped state is a bit slower, 1 s 
is still sufficient to transfer about 84 % to the strongly trapped states. This proves that after 1 s of accumulation, 
practically all molecules in the states |0;3,3,M) and |0;4,3,M) are in the top M-sublevels. 

Finally, we explain the extremely steep initial increase of the signal. This is due to the strong MW and IR fields 
coupling the three moderately trapped states. As a result, even for very short accumulation duration, one third of 
the molecules initially in the v = moderately trapped states are transferred to the state |1; 3, 3, 3). These molecules 
then have four seconds to decay back to the v = states, with a large fraction ending up in the strongly trapped 
states. 
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FIG. 8: State distribution during cooling. We show the fraction of molecules in the strongly trapped states remaining at 
the initial energy, the fraction of molecules in the moderately and weakly trapped states whose energy has been reduced by /rf 
and 2/rf, respectively, and the fraction of molecules which have experienced at least one spontaneous decay to a more strongly 
trapped set of states, and have therefore been cooled. Note that the populations in the moderately and weakly trapped states 
have been scaled by a factor of 3 and 10, respectively. 



2. Cooling 

We now apply the rate-equation model to opto-electrical cooling. This is considerably more difficult since the 
energy of the molecules now plays a substantial role. In particular, the RF transitions not only change the energy of 
the molecules but also occur at a rate strongly dependent on the energy, forcing us to keep track of the molecules' 
energy in addition to their internal state. One approach to take into account the energy of the molecules is to switch 
from discrete state populations pi to energy distributions for each of the molecular states Pi(E). This approach 
can be simplified considerably by making the following approximation. Referring to appendix [Bj we approximate 
the RF rate to be a constant times the relevant Clebsch-Gordan coefficient for E > M/rf and zero otherwise. 
This allows the energy of the molecules to be discretized: is the population of molecules in state i with energy 
^/rf < E < (n + 1)/rf- RF transitions only couple populations pi y7l with n > Mi where Mi is the M quantum 
number of the state i. Moreover, each RF transition changes n by 1. 

Finding a suitable cooling measurement for comparison to our model is considerably more difficult than for accu- 
mulation due to the inability to measure state-dependent energy distributions of the molecules sufficiently accurately. 
We therefore perform an indirect measurement, the number of molecules cooled as a function of the RF power. Ex- 
perimentally, we perform the full cooling sequence to 390 MHz. However, the RF power of the final RF at 390 MHz is 
varied. This is followed by a 1 s RF knife-edge filter at 390 MHz to eliminate molecules which have not been cooled. 
The detected signal as a function of the final RF power for cooling is shown in Fig. [7[ 

As discussed in appendix |Bj each 4 s interval in the cooling sequence is expected to result in approximately one 
cooling cycle. As a result, we chose the initial state distribution for the model to consist of molecules in the strongly 
trapped states with 2 /r F < E < 4/r F - The population with E < 3/rf is unaffected by the RF and provides a 
constant background. The rest of the population can be cooled by one additional cycle. For comparison with the 
experimental data, we show the calculated number of molecules with E < M /rf after 4 s of cooling as a function of 
the RF power in Fig. [7| Such molecules are no longer affected by the RF and would therefore also be unaffected by 
the knife-edge filter. As fit parameters we have the initial number of molecules with E < 3 /rf and E > 3 /rf as well 
as a scaling factor for the RF power. 

Although the agreement between the model and the experimental data is by no means perfect, the model and the 
experimental data show the same key features. In particular, the signal strongly decreases for both too much as well 
as too little RF power. The decrease in signal for low power is a key signature of our cooling process: as the RF power 
is ramped up, RF transitions to lower M-sublevels start to occur, but the IR rapidly pumps these molecules back to 
the top M-sublevels leading to cooling. As a result, the molecules no longer have enough energy to be pumped to 
untrapped M = states by the much stronger RF knife after the cooling sequence, resulting in a signal increase. The 
large factor 2 signal difference between optimal RF power and no/too much RF power is remarkable, indicating the 
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large number of molecules which perform at least one cooling cycle. Note that similar measurements for all other RF 
frequencies in the cooling sequence result in a similar signal increase, and have in fact been used to determine the 
optimal RF power for all RF frequencies. 

As a final application of the rate-equation model to cooling, we show the calculated time evolution of the population 
of molecules in the various molecular states in Fig. [8] We again start with all molecules in the strongly trapped states, 
but with energy 3/rf < E < 4/rf, allowing one final cooling cycle. Note that for all cooling measurements where 
constant RF cooling frequencies are applied for 4 s at a time, we increase the RF power by 2 dB after 2 s to increase 
the cooling efficiency. Similarly, we here increase the RF rate by 2 dB after 2 s resulting in the kinks at t = 2 s. The 
initial RF rate is chosen such that the largest number of molecules has been cooled after 4 s. 

The main result of Fig. [8] is to confirm that throughout cooling, almost all the molecules remain in the strongly 
trapped states, with the relative population in other states never exceeding 15 % and with the population substantially 
lower by the end of the 4 s of cooling. Even taking into account the simple nature of our model, we can therefore 
expect the great majority of molecules in the states |0;3,3,M) and |0;4,3,M) to end up in the top M-sublevels for 
all measurements shown in the main paper except no-cooling. 

Despite the small fraction of molecules in the less strongly trapped states, the calculations for Fig. [8] still show that 
approximately 9% of molecules are lost to untrapped states. We compare this value to the losses observed during 
cooling. A 9% reduction in molecules every 4s results in a factor 1.46 fewer molecules during the final 16s of the 
cooling cycle. Together with the losses due to the finite trap lifetime, these losses account for a factor 2.8 fewer 
molecules for the cooling cycle to 390 MHz compared to the cooling cycle to 1220 MHz, which is practically the same 
as the losses observed in the experiment. The losses predicted by the rate-equation model are thus consistent with 
the losses observed during opto-electrical cooling. 



Appendix D: Internal-State-Discriminating Measurement 

In the experiment, molecules are detected with a quadrupole mass spectrometer which is insensitive to the internal 
molecular state. We therefore perform microwave depletion (MWD) measurements to eliminate molecules in the 
rotational states involved in the cooling, allowing the remaining background from molecules in other rotational states 
to be subtracted. After a detailed description of the MWD, we show the results from a MWD saturation measurement 
to prove the effectiveness of the MWD. The nonzero saturation level suggests the existence of decay channels from 
the vibrationally excited state to states other than J = 3 and J = 4. We also discuss the contribution of molecules 
with | if | = 1 and \K\ = 2 to the MWD signal. 

For MWD, we apply a MW pulse to the molecules to drive all AM = ±1 transitions between the states |0; 3, 3, M) 
and |0; 4, 3, M) with M > 0, as shown in Fig. [9^. This mixes the populations in all these states. Molecules which reach 
the M = states are lost from the trap practically instantaneously. The effect of the MWD on the trap unloading 
signal is shown in Fig. |9Jd. The difference in signal with and without the MWD constitutes the state-discriminated 
unloading signal, used for the measurements in the main paper. 

Fig. |9j3 shows the number of molecules surviving the MWD pulse as a function of the duration of the pulse. A 
MWD pulse duration of 4 s, used for the MWD measurements in the main paper, is clearly sufficient to obtain 
strong saturation of the depletion. Since the relative transition rates for all the MW transitions are known from the 
Clebsch-Gordan coefficients, the differential Stark shifts, and the relative MW intensities, the rate-equation model 
from appendix [C] should accurately describe the data, predicting the solid curve in Fig. [9)3. The model allows us to 
extract absolute transition rates for all microwaves applied in the experiment. 

In principle, MWD can also be used to resolve the M-sublevels of the rotational states by driving only a subset of 
the transitions shown in Fig. [9^i. However, attempts to be M-resolved suffer from Stark broadening which prevents 
perfect separation of the MW transitions. In particular, the |0; 3, 3, M) to |0; 4, 3, M — 1) transitions have to be driven 
with high MW power due to both their small Clebsch-Gordan coefficients and their high differential Stark shifts. Each 
of these transitions therefore inevitably also drives other MW transitions in electric fields outside the trap centre, 
precluding simple M-resolved measurement. However, the results from appendix [C] prove that after both cooling and 
accumulation, the molecules are predominantly in the strongly trapped states |0;3,3,3) and |0;4,3,4). 

Despite clear saturation of the MWD, a substantial fraction of the molecules remain unaffected by the MWD even 
after 8 s. This is surprising considering the laser-off signal in Fig. [9}d is substantially lower than the signal with 
the MWD and the laser applied: the remaining molecules must have been addressed by the laser at some point. 
Considering that the |0;3,3,M) to |1;3,3,M) transition is relatively isolated in the v\ vibrational band, the most 
plausible explanation for the remaining molecules is that they have decayed to an excited vibrational ^5 state. This 
is possible due to the strong Fermi-resonance between v\ and 2^5, as discussed at the end of appendix [Aj However, 
precisely determining which states the remaining molecules occupy would require new MW or IR sources. 

Although the MWD pulse is optimized for \K\ = 3, the zero-field frequencies of the |0;3, 1,M) to |0;4, 1,M) and 
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FIG. 9: Microwave Depletion, a, Set of microwave transitions applied to the molecules during the MWD pulse to remove 
molecules in the rotational states |0; 3, 3, M) to |0; 4, 3, M) from the trap. Av/vst = j{j~fi) denotes the Stark detuning of the 
rotational states in the homogeneous-field central region of the electric trap relative to the Stark shift z/st = |E| • |d e i| = 720 MHz 
of a maximally aligned molecule. Note that transitions with AM = dzl are favoured since the microwaves are expected to be 
polarized perpendicular to the electric field in the center of the trap, b, Unloading signal after the complete cooling cycle as 
well as for the same cycle but with either the IR laser left off or the MWD applied for 4 s. The shaded region represents the 
state-discriminated unloading signal, c, Variation of the remaining molecule signal as a function of the MWD pulse duration. 
The MWD pulse is applied during a fixed 8 s of storage after the complete cooling cycle. The solid curve shows the predicted 
dependence based on rate equations. 



|0; 3, 2, M) to |0; 4, 2, M) transitions are less than 30 MHz detuned from the zero-field frequency of the |0; 3, 3, M) to 
|0;4, 3,M) transition. As a result, the MWD pulse also drives these transitions, and molecules in the states J = 3 
and J = 4 with \K\ = 1 and \K\ = 2 will also contribute to the state-discriminated signal. This effect will be largest 
for the no-cooling reference measurement since the population of molecules in the states with \K\ = 3 has not yet 
been enhanced by accumulation or cooling. Nonetheless, the fraction of molecules loaded into the trap in the states 
with J = 3, 4 which are in the states \K\ = 1,2 is only about 25 % due to reduced Stark shifts and half the hyperfine 
degeneracy. The density measured for the no-cooling is therefore overestimated by only a small amount. 

For the accumulation and cooling measurements, the contribution to the state-discriminated signal of molecules 
with \K\ = 1, 2 is considerably reduced. In particular, the zero- field A J = IR transitions from the states |0; 3, 1, M), 
|0;4, 1,M), |0;3,2,M), and |0;4,2,M) are all 4 GHz or more blue detuned from the |0;3,3,M) to |1;3,3,M) transi- 
tion [30], and the laser is therefore red detuned by a similar amount. As a result, the only IR transitions from these 
states which can occur, if they occur at all, are AM = — 1 transitions, resulting in optical pumping to untrapped 
M = states. From this we see that, compared to no-cooling, the absolute number of \K\ = 1,2 molecules in the 
accumulation signal will, if anything, decrease and the relative number is reduced by at least a factor of 2 due to the 
increased number of molecules with \K\ =3. For the cooling measurements, the RF fields will lead to strong losses for 
molecules with \K\ = 1, 2, strongly reducing their absolute number. In fact, for the complete cooling cycle to 390 MHz, 
the molecules with \K\ = 1,2 must be entirely contained in the "laser-off" signal, and since this is practically zero, 
the state-discriminated signal essentially only includes molecules with \K\ =3. 
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Appendix E: Details of the Experiment 

IR radiation at 2966 cm -1 to drive the J = 3, A J = transition of the v\ vibrational band is produced by a CW 
optical parametric oscillator locked to a frequency comb. Line assignments and frequency values published in Ref. [30] 
are verified with saturation spectroscopy with sub-MHz accuracy using a multi-pass Herriot cell and validated via 
combination differences between the \v\ J, K) and \v\ J + 1, K) states. The 410 mW of IR power illuminates the trap 
from the side (as shown in Fig.[l]3 in the main text), driving the AM = +1, A J = transitions shown in Fig. [5] with 
an estimated rate on the order of 100 Hz. 

To couple the v = 0, J = 3 and J = 4 states, MW radiation at 204.24 GHz is generated by duodectupling (xl2) the 
output of a frequency synthesizer. On the order of 5 mW MW power radiate from a MW horn antenna onto the front 
of the trap. During the cooling cycle the effective power is reduced by applying the MW with an overall duty cycle of 
25%, providing sufficient MW intensity to alternatingly drive the three narrowband J = 3 to J = 4 transitions with 
an approximate rate of 100 Hz. For efficient MW depletion we require the full 100% duty cycle. 

RF to couple neighbouring rotational M-sublevels and to apply the RF knife filter is applied directly to the 
microstructure electrode contact leads. Naturally occurring electric resonances are exploited to inject sufficient RF 
power. 

Due to the varying electric fields in an electric trap, transitions between various M-sublevels will be addressed by 
the same MW or IR frequency at different locations in the trap, potentially driving unwanted transitions to untrapped 
states. The choice of transitions shown in Fig. [5] eliminates this problem for the IR and strongly reduces it for the 
MW. In particular, for the IR transitions, the M-sublevels of a given rotational state have very similar Stark shifts 
for both the v = and the v = 1 state. In non-zero fields, all AM = — 1 transitions occur at a red-detuned frequency 
compared to the zero- field frequency, the AM = transitions occur near the zero- field frequency, and the AM = +1 
transitions are blue detuned. We can therefore easily avoid the AM = — 1 and AM = transitions while driving all 
transitions between the states |0; 3, 3, M) and |1; 3, 3, M + 1) simultaneously by blue detuning the IR laser from the 
zero-field frequency. 

Driving the three MW transitions from |0;3,3,M) to |0;4, 3, M + 1) shown in Fig. 5] while avoiding unwanted 
MW transitions turns out to be particularly easy. Due to the small differential Stark shifts between these pairs of 
states as well as the relatively large Clebsch-Gordan coefficients, relatively low MW power is sufficient to drive these 
transitions. Due to the substantially larger differential Stark shifts of the unwanted transitions, the low MW power 
and the homogeneous electric fields in the trap cause the unwanted transitions to occur at a sufficiently slow rate. 

Due to low detection efficiency, extracting a velocity distribution using a time-of-flight measurement requires an 
optimized data acquisition strategy to increase the signal per cycle. Therefore, instead of using only the first rising 
edge of the unloading signal, the quadrupole guide to the QMS is repeatedly switched on and off during each trap 
unloading. This allows several time-of-flight profiles (rising edges) to be recorded each experimental cycle which are 
summed up. In addition to improving the data acquisition rate, this results in a rising edge which represents an 
average over the velocity distribution of molecules leaving the trap at different times after the unloading process has 
started. 

To determine the density of the molecular ensemble in the trap, a density calibration measurement has been 
performed. This calibration is necessary since the molecular cloud is not measured directly in the trap but transferred 
to the ionization volume of the QMS via a quadrupole guide, which samples the molecule density in the trap. The 
flux of molecules from the quadrupole guide is obtained by integrating the density of molecules measured by the QMS 
over the transverse position of the QMS and multiplying with the velocity of the molecules. A lower bound for the 
number of molecules in the trap is then obtained by integrating the flux over time during unloading. This value 
assumes that all molecule losses occur via the exit guide. Accurate relative values for the density in the trap can be 
obtained by comparing the initial molecule density at the QMS with the size of the exit hole from the trap, which is 
proportional to E 2^ os0 ^2 • T is the temperature of the molecules, E is the unloading electric field strength, and (cosO) 
is the state-dependent orientation of the molecules. 



Appendix F: Unloading Electric Field Strength 

As already discussed in the main text, the unloading electric field strength plays a key role in properly measuring 
the velocity distribution of the molecules in the trap. Specifically, slower molecules are not efficiently extracted from 
the electric trap at high unloading fields whereas faster molecules are lost from the experiment at lower unloading 
fields due to weaker confinement. To investigate this effect, we have measured the velocity distribution for cooling 



to 1.22 GHz for a variety of unloading fields, with three representative distributions shown in Fig. 10 As expected, 



an unloading field which is too low results in a strong decrease in fast molecules, and an unloading field which is 
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FIG. 10: Effect of the unloading electric fields on the measured velocity distribution. We show the velocity 
distribution for cooling to 1.22 GHz (16kV/cm) from the main paper as well as the distribution for both half (8kV/cm) and 
twice (32kV/cm) the unloading electric field strength. The correct unloading electric field results in a measured velocity 
distribution representing practically all molecules in the trap. 

too high reduces the number of slow molecules. The unloading field of 16kV/cm with the maximal signal, chosen 
for the main paper, results in a velocity distribution which encompasses the other two distributions almost entirely, 
demonstrating that this field strength has been suitably chosen. For cooling to 390 MHz, the result of varying the 
unloading field about the optimum is practically the same as for 1.22 GHz. For accumulation, the increase in slower 
molecules for lower fields and the increase in faster molecules for higher fields is slightly larger, with the behaviour 
otherwise similar, indicating a slightly wider velocity distribution in the trap which is difficult to measure with a 
single unloading voltage. We emphasize that in all three cases, the unloading field has been varied about the optimal 
values, which vary by over a factor of 10, and not about a fixed value. The present measurement thus confirms the 
necessity to adjust the unloading fields according to the temperature of the molecules. 



